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Abstract 

We derive an analytic expression for one-loop effective action of QCD-I-QED at zero and finite 
temperatures by using the Schwinger proper time method. The result is a nonlinear effective action 
not only for electromagnetic and chromo-electromagnetic fields but also for the Polyakov loop, and 
thus reproduces the Euler-Heisenberg action in QED, QCD, and QED-I-QCD, and also the Weiss 
potential for the Polyakov loop at finite temperature. As applications of this “Euler-Heisenberg- 
Weiss” action in QCD-I-QED, we investigate quark pair productions induced by QCD-I-QED fields 
at zero temperature and the Polyakov loop in the presence of strong electromagnetic fields. Quark 
one-loop contribution to the effective potential of the Polyakov loop explicitly breaks the center 
symmetry, and is found to be enhanced by the magnetic field, which is consistent with the inverse 
magnetic catalysis observed in lattice QCD simulation. 
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I. INTRODUCTION 


The very first stage in a high-energy heavy-ion collision is dominated by extremely strong 
chromo-electromagnetic (chromo-EM) helds reflecting colliding nuclei hlled with high-density 
gluons (color glass condensate). Such a state with strong helds is called a “glasma” which is 
named since it is a transitional state between a color glass condensate (before the collision) 
and a quark-gluon plasma (QGP) [T]. The glasma is characterized by a held strength T of 
the order of the saturation scale: gT ~ Q\ (with g being the QCD coupling). Notice that 
the saturation scale is a semihard scale representing a typical transverse momentum of 
gluons in a colliding nucleus and can become large enough, at high energies, compared to 
light quark masses Qs 3> mq. Besides, it has long been known that heavy-ion collisions, 
with electrically charged nuclei, are accompanied by electromagnetic (EM) helds, but only 
recently was it seriously recognized that the strong EM helds could ahect time evolution 
of heavy-ion collision events since the strength F of the EM helds could be as large as 
or even greater than the nonperturbative QCD scale Aqcd, namely eF ^ Aq^j^ and thus 
eF S> Since both the chromo-EM and EM helds created in heavy-ion collisions 

can be strong enough compared with the light quark masses, the ehects of strong helds 
cannot be treated as perturbation (even though the coupling constants are small), but 
must be treated in a nonperturbative way. Then we expect nonlinear and nonperturbative 
phenomena associated with the strong helds to occur. Typical examples of such phenomena 
include particle productions (quarks, antiquarks and gluons) from these strong helds (the 
Schwinger mechanism), which must be a key towards understanding the formation of QGP. 

While the (coherent) chromo-EM helds will disappear as the QGP is formed, the EM helds 
could survive longer due to Faraday’s law, which works in the presence of a conducting 
medium mm- If the EM helds survive at a strong enough level until the formation of 
QGP, and even until the end of the QGP’s lifetime, we need to describe the QGD phase 
transition with the ehects of strong EM helds taken into account. Notice that the ehects 
of strong magnetic helds on thermodynamical or fundamental quantities of QGP can be 
investigated in lattice QGD simulations, and are indeed found to be large. For example, at 
zero temperature, lattice QGD simulations conhrmed the “magnetic catalysis” as predicted 
in several ehective models [81-IT5] in which the value of chiral condensate increases with 
increasing magnetic held strength. On the other hand, at hnite temperature, lattice QGD 
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simulations almost at the physical point concluded [161 dZ] that the magnetic catalysis 
does not necessarily occur at all the temperature regions, but rather gets weakened and 
even shows opposite behavior with increasing temperature. Such behavior of the chiral 
condensate around the critical temperature is called “magnetic inhibition” [18] or “inverse 
magnetic catalysis”, which eventually gives rise to decreasing critical temperature. For recent 
reviews on the phase diagram of chiral phase transitions in strong magnetic helds, see, e.g.. 
Refs. [201 121]. Furthermore, it is reported [19] that the (pseudo)critical temperature of 
the conhnement-deconhnement phase transition (for the Polyakov loop) also decreases with 
increasing magnetic held. This is achieved by increasing Polyakov loop expectation values. 
Probably, these two phenomena are related to each other. However, so far, there is no clear 
explanation about the physical mechanism behind this (for recent attempts, see Refs. [22l 123] 

and mm)- 

We can investigate these two aspects, namely the nonlinear and nonperturbative dynam¬ 
ics of strong helds (including particle production) and the phase transition under strong 
external helds, within a single framework of an ehective action. So far, ehective actions 
for QED and QCD in various external conditions have been extensively explored. First of 
all, Euler and Heisenberg derived a nonlinear ehective action for constant EM helds at the 
electron’s one-loop level, known as the Euler-Heisenberg (EH) action [26]. Later, Schwinger 
reproduced the same action in a held-theoretical manner, which is the so-called Schwinger 
proper time method The EH action at hnite temperature is computed in imaginary 
time formalism [2S1122] as well as in real time formalism [301 El]- Furthermore, an analog 
of the EH action in QCD (for chromo-EM helds) has been evaluated too within a similar 
method at zero and hnite temperatures mm- Lastly, the most recent progress was to 
compute the EH action at zero temperature when both the EM and chromo-EM helds are 
present, which was done by one of the authors and B. V. Galilo and S. N. Nedelko in¬ 
dependently [321113] • The author of Ref. [13] used this ehective action to investigate the 
QCD vacuum (gluon condensate) in the presence of strong magnetic helds. Though all of 
these are about the ehective action for strong helds and choromo-EM condensates, it should 
be possible to include the Polyakov loop at hnite temperature. Indeed, an ehective action 
(or potential) for the Polyakov loop at the one-loop level was computed independently by 
D. J. Gross, R. D. Pisarski, and L. G. Yahe [31], and by N. Weiss [331133], and the result 
is called the Weiss potential. In the present paper, we are going to derive an analog of the 
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EH effective action in QCD+QED at finite temperature with the Polyakov loops included. 
Thus, the result may be collectively called the “Euler-Heisenberg-Weiss action.” Our result 
is also a generalization of the one obtained by H. Gies [H], who computed an effective action 
for the Polyakov loop and the chromo-electric held. 

The paper is organized as follows: In the next section, we will derive the effective action 
for QCD-I-QED at hnite temperature by using the Schwinger proper time method. Variables 
of the effective action are the EM and chromo-EM helds as well as the Polyakov loop, and 
one can reproduce the previous results (the EH action with QCD-I-QED helds, the Weiss 
potential, etc.) in various limits. Then, we discuss some applications of our ehective action 
in Sec. HI. First, we investigate quark-antiquark pair production in QCD-I-QED helds at 
zero temperature. We obtain the quark production rate in the presence of QCD-I-QED 
helds, which allows us to study the quark pair production with arbitrary angle between the 
EM and chromo-EM helds. Next, we study an ehective potential for the Polyakov loop 
with electromagnetic helds. We hnd that the magnetic held enhances the explicit center 
symmetry breaking, while the electric held reduces it. This indicates that the (pseudo)critical 
temperature of the conhnement-deconhnement phase transition decreases (increases) with 
increasing magnetic (electric) held. Finally, we conclude our study in Sec. IV. 

II. ONE-LOOP EFFECTIVE ACTION FOR QCD-hQED AT FINITE TEMPERA¬ 
TURE 

In this section, we derive the one-loop ehective action for QCD-I-QED at hnite temper¬ 
ature. The ehective action will be a function of chromo-EM and EM helds, as well as the 
Polyakov loop. Notice that both the strong helds and the Polyakov loop can be treated as 
background fields so that the background held method is applicable. We will take quantum 
huctuations around the background helds up to the second order in the action, and integrate 
them in the path integral. This corresponds to computing the action at the one-loop level. 

We shall begin with the four-dimensional QCD action of the SU{Nc) gauge group with 
Nf havor quarks interacting with EM helds: 

Sqcd+qed = + (1) 
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where the covariant derivative contains gluon helds^ ^4“ (a = 1,..., — 1) and U(l) gauge 

helds as 

D^ = d^- - ieQqttf ,, ( 2 ) 

and the gluon and EM held-strength tensors are given by + gf°‘^'^A^^A’l 

and = d^Qy — dya^ , respectively. In this paper, we treat the EM fields just as background 
fields, and assume that the field strengths are constant so that Of = 0. We abbreviate color, 
fiavor, and spinor indices of the quark field in Eq. Q. Mass and charge matrices of quarks 
are given by Mq = diag(mq^, ..., niq^^) and Qq = diag(Qq,, ..., Qq^,^). As for the 
gluon field, we apply the background field method and decompose the gluon field into a 
slowly varying background field and a quantum fluctuation as 

Al = Al + Al. (3) 

Here we employ the covariantly constant field as a background field, which obeys the fol¬ 
lowing condition 


= 0 . ( 4 ) 

where the covariant derivative is defined only with respect to the gluon background field: 

= ( 5 ) 

and = d^Ay — dyA^^ +gf‘^^'^A!‘^Al. From the condition (|^, the field-strength tensor 
can be factorized as where rA is a unit vector in color space, normalized as 

= 1, whereas Tf^y expresses the magnitude of the chromo-EM field. We further assume 
that J^fj,y is very slowly varying, satisfying d^J^^y = 0, which allows us to obtain the analytic 
expression of the EH action for QCD, just as in QED. Both T^y and n“ are space-time 
independent. The background field is proportional to the color unit vector n“ as 

= AX, (6) 

^ Throughout the paper, we use a, 6, c (and h) for adjoint color indices {a,b,c = — 1), i for 

fundamental color indices {i = l,...,iVc), for Lorentz indices, and / for flavor indices (/ = 

1,...,W)- 
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and the field-strength tensor has an Abelian form, = d^^Ay — dyA^. This backgronnd 
field (|^ indeed satisfies the condition (|^. By nsing the backgronnd field and the qnantnm 
finctnation, the fnll glnon field-strength tensor can be decomposed as 

= F,yn^ + - iP-ip + gr^^A^Al. (7) 

Applying the backgronnd gange for the qnantnm finctnation, 


v^:ai = 0, 


( 8 ) 


we get the gange fixed action in the presence of EM fields. 


*S'qcd+qed 



-\+ (yrT - ®rT)+ 9 /“‘AT}" 

-? {V^D^r c‘ + q dl^D^ -M,)q- , 


(9) 


where c is the ghost field and ^ is the gange parameter. Notice that one of the covariant 
derivatives in the ghost kinetic term D^f and the one in the qnark kinetic term 77^ defined 
in Eq. ([^ contain all the gange fields. The effective action for the backgronnd fields and 
can be obtained throngh the fnnctional integral as 


exp( zS'efr[M;_j, a^]) = ^A^c^cS>qS>q exp (z 


d a;S'QCD+QED 1 


( 10 ) 


We perform the fnnctional integral with finctnations taken np to the second order. This 
corresponds to evalnating the one-loop diagrams as shown in Fig. 1. The glnon, ghost, and 
qnark loop integrations can be separately done, and one finds, respectively. 


^iexp / - ‘lgr^^Fly\ A^'^ \ = det[-(p2)-^^^ _ - , 


( 11 ) 


+1 


exp <j 7 / dS c“ [-{Vy^] A \ = det , 


^q^q exp 7 z / d'^x q (zq^T^^ — Mq ] q \ = det — M, 


Here we have taken the Feynman gange, .^ = 1. In the qnark one-loop contribntion, the 
covariant derivative contains both of the backgronnd fields A^ and a^\ 


F/j, Ffi icQqClqi 

= dq,- igAlT'^ - ieQqa^ . 


( 12 ) 
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gluon loop 


ghost loop 


quark loop 


FIG. 1: Typical loop diagrams contributing to the effective action. The field A contains both the 
chromo-EM fields and the Polyakov loop. 

On the other hand, the gluon and ghost one-loop contributions contain and which 
only depend on the gluon background held A^. This is, of course, because the gluon and 
ghost helds do not have electric charge and thus cannot interact with EM helds. Since these 
contributions are the same as in the pure Yang-Mills (YM) theory, we may call these the 
YM part. 

So far, we have not specihed the background held but it can contain both the chromo- 
EM helds and the Polyakov loop. Let us briehy explain how the Polyakov loop is described 
within our framework. In the pure Yang-Mills theory at hnite temperature, there is a 
conhnement-deconhnement transition whose order parameter is given by the Polyakov loop. 
It is dehned by the (closed) Wilson line along the imaginary time (r) direction: 

= ^TrPexpji^^ dr/l 4 (r,f)T“| , (13) 

where /? = 1/T is the inverse temperature and V stands for a path-ordered product along the 
imaginary time direction. Indeed, (<I>) —)■ 0 (($) 7 ^ 0) corresponds to a conhning (deconhned) 
phase, since the negative logarithm of the expectation value of the Polyakov loop can be 
identihed with the free energy of a static quark (a vanishing value of the Polyakov loop 


7 








implies that the energy of a single quark state is inhnity). These two phases are distinguished 
by the center symmetry. The gauge helds at hnite temperature are not necessarily periodic 
in the direction of imaginary time and can have ambiguity related to the center subgroup 
of the gauge symmetry SU(A'"c). This residual symmetry is called the center symmetry 
and the theory is invariant under gauge transformations which differ at r = 0 and t = (5 
by a center element of the gauge group. The Polyakov loop $ transforms as $ —)■ 

(n = 0,1, 2,..., ~ !)■ Thus, the values of distinguish the center symmetric (conhning) 

phase and the center broken (deconhned) phase. Dynamical quarks, however, explicitly 
break the center symmetry. Therefore, in QCD, the Polyakov loop should be understood 
as an approximated order parameter. Still, we can compute an effective action for the 
Polyakov loop and discuss how a phase transition occurs when external parameters such as 
temperature are varied. 

An effective action for the Polyakov loop in the pure Yang-Mills theory was obtained in 
Refs, im 03] in the following way: Working in what we now call the “Polyakov gauge” for 
a time-independent held Al{x) = in the SU(2) case, the authors of Refs. [lH US] 

performed a functional integral with respect to huctuations around the held 0(5;). This 
procedure is nothing but the one we explained above where we treated the gluon held as 
a background with a huctuation around it. Besides, as long as we consider a spatially 
homogeneous and time-independent order parameter we can have both the Polyakov 
loop and the chromo-EM helds at the same time. We divide the background held into 
the constant part and the coordinate-dependent part as + A^{x))rf‘. The 

second term gives the real (physical) chromo-EM helds so that = dfj_AAx) — dyA^^{x) = 
(0^Mi.(a;) —0,yM^(x))n“, while the hrst constant term does not. We want to treat both the 
chromo-EM helds and the Polyakov loop, and the latter is described at hnite temperature. 
In order to have the both, we specify the transformation of the temporal component of the 
background held Aq{x) under the Wick rotation of the coordinate, Xq —)■ —ix^ = —ir and 
Xi Xi {i = 1 , 2, 3), as follows: Aq{x) = (Mo + Mo(a;))n“ —)■ {iAi + Ao{x))n°'. In this way. 


the hrst term gives the Polyakov loop dehned in Eq. (13), while the second term remains 



unchanged to give the real chromo-EM helds. We work in the Polyakov gauge for A, 


a rr^2. 


^^4 — ^^4 6 ^— 0 , 


(14) 


which does not conflict with the covariantly constant condition in Eq. (|^. Notice that we 
use this gauge with 5“^ even for the SU(A^c) case, and the color unit vector introduced 
in Eq. ([^ should be understood as = 6^°’ at hnite temperature.^ Following Ref. |13], we 
also introduce a dimensionless held C as 


^ _ 9^4 
27rT’ 


so that the Polyakov loop is simply given as 


(15) 


$ = cos(7rC) for SU(2), 

$ = -|l + 2cos(7rC')| for SU(3). 


(16) 


A. Yang-Mills part of effective action 


Now, we consider the Yang-Mills part (gluon and ghost contributions) of the one-loop 
effective action. In the one-loop level, the effect of EM helds is not included in gluon and 
ghost loops, since these do not directly interact with EM helds. From Eq. 0 , the ehective 
actions of gluon and ghost parts are given, respectively, as 


= In det [-{vy^g,, - ^ , 

(17) 

z^ghost = In det [-( 1 ) 2 )-]+i_ 

(18) 


^ In the literature, the fourth component of the gauge field A 4 in the Polyakov gauge is often expressed 
in terms of W — 1 real scalar fields. In our formalism, these fields are properly encoded in the color 
eigenvalues (i = 1,..., W) and Vh {h = 1,..., — 1), which will be defined later. Here, choosing the 

third direction of the color unit vector—= 5“^ at finite temperature—we pick up the one particular 
field Ai which provides a simple expression for the Poyakov loop as shown in Eq. (16). However, in the 
finial expression of our effective action, it is quite straightforward to keep all the W — 1 scalar fields in 
the color eigenvalues coi and Vh- 

^ Still, we keep the expression n“ because we will discuss the case at zero temperature. 
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Let us first explore the gluon part (17). By using the proper time integral, the gluon part 
of the effective action can be rewritten in the following form (the limit e, 5 —)■ 0 is always 
implicit and should be taken after the calculation): 




Af2-1 

h=i 

N^-l 

- K t . 

h=l 


el-e 'I I ' 


ds 


-5s ^^-i{2gvha)s _|_ ^-i{-2gvha)s j_ ^-i{igvhb)s j_ ^-i{-2igvhb)s 


+ e“ 


+ e 


x). 


(19) 


While the capital trace “Tr” in the hrst line is taken with respect to colors, Lorentz indices, 
and coordinates, “tr” in the second line is only for Lorentz indices. Also, in the second line, 
we have introduced real quantities Uh (h = 1,..., N'^ — 1) that are eigenvalues of a Hermitian 
matrix (i.e., = Vh(p°'), and Lorentz-invariant quantities a, b dehned by 


<^=\'j\jT' + (T'TY + T\ h = \\l’jT‘ + {T-fY-T\ ( 20 ) 

with the dual held-strength tensor (or equivalently, by | and 

ab = ■ J^). The covariant derivative is dehned as — igvhA^. The calculation 

up to now is in fact the same as in the case at zero temperature which was done in Ref. |13] . 
At hnite temperature, however, one needs to be careful in evaluating the matrix element 
Namely, it can be now written as the Matsubara summation: 


(a;|e 


\x) = 


n 




(27r) 


( 21 ) 


pQ=igVhA4—i27rnT 


where the functions X^^{s) and W(s) have been dehned as m 


= [(d'^hA') Han(gVhXs)]°‘^, 

^h(s) = ^tr In cos(gVkXs). (22) 


We use the following identity: 


1„(A1 - li) = ) - Ag 


ds 




in the limit e —)■ 0 and 5 —)■ 0. We ignore the first divergent term, since it does not depend on the fields. 
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In the presence of the Polyakov loop A 4 , the periodic boundary condition of the gluon in 
the imaginary time direction is modihed. Then, the Matsubara frequency is shifted by the 


Polyakov loop as in Eq. (21). Performing the three-dimensional momentum integral and 


applying the Poisson resummation [50] , one can obtain the matrix element in terms of a and 


b as 

where 

[)(s) = 

with 


gvhus 


gvhbs 


lOvr^ sin(5fr;/ias) sinh(5fn/ibs) 


00 

1-1-2 ^ e^^’^^cos I —n 

n=l ^ 


/ gvhA4 


,(23) 


b^ - 

-I- b^ 


gvhacot(gvhas) + 


-h 

-I- b^ 


gvhbcoth(gvhbs), 


(24) 




(25) 


The vector is the heat-bath four-vector, which is (1,0, 0,0) in the rest frame of the 


heat bath. The hrst (second) term in Eq. (23) corresponds to the zero-(£nite-)temperature 
contribution. The gluon part of the effective action is then given as 

W2-l 

r . 

32vr2 I - - L c3-. 


/ c pc 


^-5s ^^-i{2gvha)s _|_ ^-i{-2gvha)s ^-i{igvhb)s _|_ ^-i{-2igvhb)s'^ 


X ■ 


gvhas 


gVhbs 


sm{gvhas) smh.{gvhbs) 
Similarly, we obtain the ghost part as 


l + 2^e‘J@”’cos 


n=l 


/ gvhAj 

V T 


n 


(26) 


iS, 


ghost 


,-i+e r ds 


327r2 


d^x 


h=l 


'0 


r,3—e 


X 


gVhUs 


gVhbs 


sm{gvhas) smh.{gvhbs) 


1 + 2 x 1 


e cos 


n=l 


f gvhAj 

V T 


n 


(27) 


In both parts, the hrst terms in the square brackets are the results at zero temperature 
and agree with the known results |13]. As discussed in detail in Ref. [13|, each term has an 
ultraviolet (UV) divergence, which, however, can be absorbed by renormalizing the coupling 
g and helds [221 E3]- On the other hand, the finite-temperature contributions do not 
have UV divergence, and thus we do not need an additional renormalization procedure for 
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the finite-temperature contributions. We regard the coupling and fields as renormalized ones 


and focus on UV-finite pieces in Eqs. (26) and (27). 


Our results (26) and (27) are effective actions for chromo-EM fields as well as the Polyakov 


loop at finite temperature. These are generalizations of the previous results in two cases. 
Indeed, if we consider the pure chiomo-electric background with a Polyakov loop {B = 0, E ^ 
0, ^0 7 ^ 0), we find a —)■ b —0 and reproduce Gies’s effective action at finite temperature 
[H]. Moreover, in the case of the pure chromo-ma^'nefic background (T = 0, S 7 ^ 0, M 4 = 0), 
we find a i3, b —j- 0 and reproduce the results obtained in Refs. 


B. Quark part of effective action 

For the quark part of the effective action, we follow basically the same procedures as 
in the Yang-Mills part. From the functional integral 0 . the quark part of the one-loop 
effective action reads 


^*S*quark 111 det 




Utilizing the proper time integral, we evaluate the effective action as 


iS, 


quark 


= Tr In 


- Mq 


i=l f=l >^0 


(28) 


• • e ^f ^00 1 


where f with the field f being a linear combination of the gluon field 


ij 


ij 


and the photon field as 


= guiA^ + eQq^a^. 




(30) 


This covariant derivative ©('j can be obtained from defined in Eq. (12) with the covari- 
antly constant field employed as the background field. Here oji {i = 1,..., Nc) are eigenvalues 
of an Nc X N^. matrix n°‘T°' and satisfy® ^ — 1/2- The field-strength 

tensor can be expressed in terms of constant chromo-EM fields S, B, and EM fields E, 


® Let 0 be a diagonal matrix with eigenvalues oji, i.e., P = diag(a;i, ... ,ujnJ = Un°‘T°‘UE Then, — 

trP = n“trT“ = 0 and = tr (r“r'’)n“n^ = 1 / 2 . 
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B as [with the notation V = (\4, Vy, Vz)] 




( 0 

-£-t. 


= 9^i 


£x 

0 


£y £z \ 

Bz -a, 


-£y -Bz 0 B^ 

\ -£z By -Bz, 0 ) 


+ gQ 


‘If 


^ 0 Ex 

Ey 

Ez 

\ 


0 

tb 

1 

Bz 

By 


1 

1 

to 

0 

Bx 



V ~Bz By 

-Bx 

0 



en by Eiaij 

and Ebij 

with 


(31) 


o 


FJ/ + (F.J • F.,/)= + F?/ . b.j 


bif — - 


F‘j + (F,,/ ■ \,r - n, ■ ( 32 ) 


The dual held-strength tensor ¥^j is dehned as ¥^j = le^''^f^¥ijafs. By using Eq. (31), 
F^j = 2(a^j — bjj) and Fjj • Fjj = Aaijbij can be expressed in terms of chromo-EM helds 
and EM helds as 


nf = 2(it/-^/)> 

= ~^£i,f ■ Bij, 


(33) 


where we have dehned the combined electromagnetic helds as £ij = guJi£+eQq^E and Bij = 
gUiB + eQq^B. Taking the trace of the matrix at hnite temperature, 

we get 




OO p 




(27r 


(34) 


PQ=igiOiA4,—i'E{2n-\-l)T 


Here, the functions X“^(s) and Yjj(s) have been dehned as [50 


= [F-;tim(F,,/S)]“y 
^ W COs(F,_ys) . 


(35) 


In the presence of the Polyakov loop M 4 , the antiperiodic boundary condition for the quark 
is also modihed. Then, the temporal component of the four-momentum vector has been 


replaced by the Polyakov loop and the Matsubara frequency for a fermion in Eq. (34). The 
third part, tr is common with the case at zero temperature and was computed in 

Ref. |13]. The result is 

tr exp ( ^(T ■ Fjjs j = 4cos(aijs)cosh(bjjs). (36) 
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Now, performing the three-dimensional momentum integral and using the Poisson resum- 


mation, we hnd from Eq. (34) 
tr(x|e-'<-®?./-^"'^-^)|a:) 


sin(aijs)sinh(bijs) 


cos(ai js)cosh(bijs) 


X 


l + 25 ;](-l)”e 




cos 


n=l 


gojiA^n 

T 


(37) 


where 


1^2 _ ^2 ^2 _ j _ ^2 

W(s) = -|^7i;t^a*,/cot(aijs) +^b*jcoth(bijs), 


</ + Kf 


with 


<j + Kj 


I, = 


(38) 


(39) 


In the heat-bath rest frame, we have = (1, 0, 0, 0) and then e^j = £fj = {gUiS + eQq^E)'^. 
Therefore, the quark part of the one-loop effective action reads 

7.^ ' I 

^*S^quark 


Stt^ 


[ *^^'(ai,/s)(t’i,/s)cot(ai,/s)coth(bijs) 

J .-I Jo * 


X 


7=1 / = 1 
oo 


l + 25 ^(-l)”e: 




COS 


n=l 


goJiA^n 

T 


(40) 


As in the YM part, the hrst (second) term corresponds to the zero-(£nite-)temperature 
contribution. The zero-temperature contribution agrees with the previous result obtained 
in Ref. [43] . 

Again, the hrst term contains UV divergences. These divergences have two origins: 
QCD and QED |33]. This is because the resummed quark one-loop diagrams contain 
contributions from the diagrams with only two EM held insertions (QED) and only two 
chromo-EM held insertions (QCD). The UV divergence coming from purely QCD dynamics 
is additive to the one which we encounter in the YM part. Then, we can absorb all the UV 
divergences by renormalizing the coupling g, e and helds A^, a^. From the renormalization 
procedure at zero temperature, we have obtained the correct beta functions of both QCD 
and QED in Ref. [13]. The sum of the three parts (26), (27), and (40) may be called the 


Euler-Heisenberg-Weiss action in QCD-I-QED at hnite temperature. This result can be 
applied to several systems where strong EM helds and chromo-EM helds coexist at zero 
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and finite temperatnres. In the next section, we will show some applications of onr effective 
actions. 


III. APPLICATIONS OF EULER-HEISENBERG-WEISS ACTION IN QCD+QED 

In this section we will discnss two applications of onr resnlts. The hrst one is the qnark 
pair prodnction in the presence of both EM and chromo-EM helds. We treat the effective 
action at zero temperatnre. The second application is to investigate the effects of EM helds 
on the effective potential for the Polyakov loop at hnite temperatnre. We will discnss the 
possible implication for the inverse magnetic catalysis. 


A. Quark pair production in QCD+QED fields 


Let us hrst discuss quark-antiquark pair production in constant QCD+QED helds as an 


application of our ehective action. For this problem, only the quark part (40) is relevant. 


In the early stage of relativistic heavy-ion collisions, extremely strong chromo-EM helds 
and EM helds could coexist. Notice that the strong electric held in addition to the strong 
magnetic held could be created on an event-by-event basis [5] . The strength of the chromo- 
EM helds is approximately of the order of the saturation scale: \gB\, \g£\ ~ Q^, whereas 
strengths of EM helds would reach the QCD nonperturbative scale \eE\, \eB\ ^QCD^ 
even exceed it. Under such strong QCD+QED helds, a number of quark-antiquark pairs 
must be created through the Schwinger mechanism. The pair-production rate per unit space- 
time volume can be obtained from the imaginary part of the quark ehective Lagrangian at 


zero-temperature. Taking the zero temperature contribution in Eq. (40), one hnds 

N, Nf 


c 


quark 


*S^quark 

/ 


1 


EE 

i=i f=i■ 


I' "^\aijs){bijs)cot{aijs)coth{bijs). (41) 

0 -S 


This is the same as the result obtained in Ref. 


The imaginary part of the ehective 
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FIG. 2; Contour on the complex s plane. The contour along the real axis is inclined by an 
infinitesimal number 5 > 0. 


Lagrangian thus reads 


3m C 


1 


Nc 


quark 


Svr^ 
1 1 


r°° ds 

VV / —e-^"sin(mj s) X (aijs)(bijs)cot(aijs)coth(bijs) 

i=l f=i -^0 s 


=1 /= 

2? ^ ^ ^ ^ 1 / c3 


ds —is(irn^ +iS) f ds —is(m^—i5) 

e f ^ e 


x(aijs)(bijs)cot(aijs)coth(bijs) 


(42) 


The integrand has inhnitely many poles along the real axis [from cot (a* js)] and along the 
imaginary axis [from coth(bijs)]. With a small positive number 5 > 0, the integral contour 
along the real axis is inclined. Closing the contour in the lower half of the s plane as depicted 
in Fig. 2 and picking up the poles lying on the imaginary axis Spoies = —inTr/bij, we hnd 

N, Nf 

^Tfl /Iquark 


i=l /=! n=l ^ ‘’J 


(43) 


By using this expression, we can investigate quark-antiquark pair productions under arbi¬ 
trary conhgurations of constant chromo-EM and EM helds. The production rate per unit 
space-time volume is given by Wgg = 23m Equark- When we take Nc = Nj = 1, Q = 1, g ^ 0, 


B —)■ 0 and replace irig —)■ mg in Eq. (43), we reproduce the well-known Schwinger formula 
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for the production rate of e+e pairs in an electric field [27] : 


We+e- = 2Sm £eh = 


(eE)‘ 


E 


n=l 


(44) 


as we expected. On the other hand, in the pure chromo-electric held case, we obtain the 
same formula for quark productions derived by G.C. Nayak isg. 

1. Quark pair production in purely electric background 


First, we shall consider quark pair production in a purely electric background with van¬ 
ishing magnetic helds: B,B ^ 0. In this case, the production rate for qq pairs of havor / 
becomes 


Nc 








1 

1 — V ——'i 




2=1 


n=l 




(45) 


ij = V 0\/ = \/ + {eQq^yE^ + 2gUieQq^SEcos9£E, with E = £ = 


where = 

\/^, and 9se being the angle between E and £. For Nc = 3, the eigenvalues Ui are given 
by oji = 1/2, U 2 = —1/2, and = 0. Recall that a factor gui plays the role of an effective 


coupling between the chromo-EM held and quarks [see Eq. (30)]. Thus, a quark (or an 
antiquark) with cja = 0 does not interact with the chromo-EM held in this representation. 
Still, since there is always a coupling with the EM helds, qq production with cua = 0 is 
possible due to electric helds, i.e., bi= 3 j = \eQqfE\ ^ 0. 

Let us see the dependences of production rates on the quark mass rriq and the angle dsE- 
We hrst consider the case with light quark masses bjj. The left panel of Fig. 3 

shows the light (up) quark production rate with rriq = 5 MeV and Qq = -1-2/3. The chromo¬ 
electric held is hxed to = 1 GeV^, which is a typical value realized in heavy-ion collisions 
at RHIG and LHG, while we take several values of strength for the E held. The production 
rate increases with increasing E held, which is an expected behavior of the usual Schwinger 
mechanism, but it does not show dependence on the angle Ose, while bjj certainly depends 
on dsE- This unexpected behavior can be understood as follows: When the quark mass is 
small enough, m? bjj, we can approximate the production rate as 




'Wqfqj ~ 




f (g£)= 




2=1 


72=1 




47r3 [ 2 


(46) 
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where C(2) = vr^/h and bj = + {eQqYE"^ + 2gUieQqSEcos9£E- Notice that the 


angle dependence in b* drops ont thanks to the relations =1/2 and = 0- 

Therefore, the prodnction rate is independent of the angle 9s e- 

We next discnss the prodnction of heavy qnark-antiqnark pairs. Since the heavy qnark 
limit jnst implies that the pair creation does not occnr, we consider the case where qnark 
masses are comparable to the backgronnd held ~ bjj. This is realized for charm qnarks if 
we again take the typical valne of the chromo-electric held gS = 1 GeV^. For rric = 1.25 GeV 
and Qq = Qcharm = +2/3, the prodnction rate of a charm qnark pair is shown in the right 
panel of Fig. 3. This time, while the prodnction rate becomes small, one can see a clear de¬ 
pendence on the angle 6 s e- Both ehects (small production rate and angle dependence) come 




from the exponential factor in Eq. (45). In particular when the electric held is parallel (or an¬ 
tiparallel) to the chromo-electric held, the production rate has a maximum. Since the expo¬ 
nential factor is very sensitive to the change of bjj, the rate is largely enhanced at 6 se = 0, vr. 
Symmetric shape of the angle dependence with respect to = 7r/2 is not so trivial. Notice 
that the ehective held strengths of the combined held at 9 se = 0 and tt are not equivalent 
for a hxed value of i; namely, it is the strongest for the parallel conhguration (for uji > 0) 
bi,chi,rra{9sE = 0) = ^J{gUiY£‘^ + (eQcharm)^-^^ + 2 gUieQE aud the weakest for the an¬ 
tiparallel conhguration b*, charm(^'^E = tt) = y^{gUiyS‘^ (eQcharm)^-^^ - 2 gUJieQchi,rm£E, 
implying that pair production is most enhanced for the parallel conhguration. This 
is true for any index of i giving a positive eigenvalue Ui > 0. However, this eigen¬ 
value appears with a partner ujj having an opposite sign ojj = —oji [for SU(3) we have 
oji = —UJ 2 = 1/2], and the antiparallel conhguration gives the strongest ehective held for 
the index j, bj^charm(^£E = tt) = bj^charm(^£’E = 0). Therefore, after summing over all the 
pairwise modes i, we obtain the angle dependence symmetric with respect to 9se — 7r/2. 


2. Quark pair production in purely chromo-EM background 


Next, we investigate quark pair production under chromo-EM helds in the absence of 
EM helds. Lorentz-invariant quantities and Fjj ■ Fjj are now explicitly given as [see 


Eq. (33)] 


= 2{gui)\B^ - 8^) , F,j ■ F,j = -4(^u;,)"FHcos Osb , 


(47) 
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FIG. 3: Quark production rate as a function of the angle which stands for O^e- The left 


panel is the light (up) quark production rate, while the right panel is the heavy (charm) quark 
production rate. The chromo-electric field is fixed a,s g£ = 1 GeV^. 


where B = Vb'^, and OgQ stands for the angle between S and B. When Osq = ±7r/2 and 
S > B, we can move into a system with pnre chromo-electric helds with = a* = 0 and 
bij = bi = Igujily/S"^ — B'^ by the Lorentz transformation. Then, the prodnction rate for a 
certain flavor of qnark becomes 


1 °° 1 rni 


29m £quark - ^^3 ^ ^ ^2 ^ 

2=1 n=l 

which decreases as B increases. Fnrthermore, for B > £ the prodnction rate vanishes since 
in this case the system is eqnivalent to the pnre chromo-magnetic held system. When 
= 0,7r, which wonld be relevant conhgnrations for relativistic heavy-ion collisions, Ui 
and bi become a* = \gUiB\, bj = \gui£\. Then, the prodnction rate reads 


] _ /jR \ 

29m £quark = ^ ^ \g(^iB\\guji£\ ^ - e'^^^'^coth i -mi j . (49) 

This prodnction rate is the same resnlt as obtained in Refs. [531154] . It increases as either the 
chromo-electric held or the chromo-magnetic held increases. Fignre 4 shows Osb dependence 
of the light qnark prodnction rate with a hxed valne of the chromo-electric held, g£ = 1 
GeV^. The maxima appear when the chromo-magnetic held is parallel (or antiparallel) to 
the chromo-electric held. 
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FIG. 4: Light (up) quark production rate as a function of which stands for with 

vanishing electromagnetic fields. We take the strength of the chromo-electric field as = 1 GeV^. 

3. Quark pair production in a glasma with EM fields 


Now we shall consider a specific configuration of chromo-EM fields that are relevant 
for relativistic heavy-ion collisions accompanied by EM helds. Suppose that the chromo¬ 
electric held and the chromo-magnetic held are parallel to each other, B \\ and that these 
strengths are approximately equal to the saturation scale: \gB\ = \g£\ = 1 GeV^ ~ 

This conhguration of chromo-EM helds is indeed realized at the very early stage of the 
glasma evolution. Under this condition, we investigate light (up) quark productions with 
niq = 0.5 MeV and Qq = -1-2/3. Let us turn on the EM helds. In the heavy-ion collisions, 
the dominant EM held is the magnetic held perpendicular to the beam direction (equivalent 
to the direction of the glasma helds). But here we consider the case \eB\ 0 and \eE\ = 0, 


with arbitrary orientation. Then, the quantities ¥jj, and Fjj • Fjj read [see Eq. (33)] 


¥lj = 2 [{eQq)^B^ + 2guieQqBBcos9BB] , 

¥ij-¥ij = -4[{gUiy£B + gUieQq£BcoseisB] , (50) 

with B = Here we have used the fact that cosO^b = cosObb- Note that in the case 

of antiparallel conhguration of B and £, results are the same as those of the parallel case, 
since this changes Fjj ■ Fjj —)■ —Fjj ■ Fjj, but it is squared in and bij. 
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Figure 5 shows the quark production rate as a function of the angle 633 
with several strengths of the magnetic held. At the angle relevant for relativis¬ 
tic heavy-ion collisions, 633 = 7 r/ 2 , the production rate slightly decreases with in¬ 
creasing B held. This can be understood from Eq. (43) as follows: In this 


case, the quantity A[eQqYB^ -|- lQ{guJiY£B + 2 {eQqYB‘^ (or bjj = 

+ lQ{guiYSB - 2 {eQ,fB‘^ ) increases (decreases) with increasing B held, 
while the product cnjbij = \£ij ■ Bij\ = {goOiYSB is independent of B held. Therefore, at 
^33 = 7 ^/ 2 , the quark production rate monotonically decreases due to the exponential factor 
exp{ —(mq/bjj)n 7 r}. This result is independent of the sign of cuj. 

On the other hand. Fig. 5 shows that the quark production rate increases with increasing 
B held at 633 = 0 and tt. This can be understood as follows: At 633 = 0 , 7 r, the quark 


production rate reads from Eq. (|43|) 

1 


Nc 


2 Am C, 


quark 


47J-2 


^\gu,\£Bij^-e 


— 0 

n 


coth 




\guJi\£ 


mi 


(51) 


2=1 72=1 

where the strength of the combined magnetic held has been dehned as Bij = \guJiB + eQqB\ 
for 633 = 0, whereas Bij = \guJiB — eQqB\ for 633 = vr. This production rate has a 


similar form with Eq. (49). First, we consider the case \goJiB\ > \eQqB\. When the chromo- 
magnetic held and the magnetic held are (anti)parallel to each other, 633 = 0 {633 = tt ), 
with oji > 0 {oji < 0 ), the strength of the combined magnetic held B^j linearly increases 
with increasing B held, and thus coth ^ |^’|g W 7 rj slightly decreases and approaches unity. 
When 633 = 0 {633 = tt) with cjj < 0 (ui > 0), the held strength Bij linearly decreases with 
increasing B held, but coth f increases. Then, after summing over all the modes 


i, the production rate (51) at 633 = 0 {633 = tt ) monotonically increases with increasing B 


held. In the case of \guJiB\ < \eQqB\, the production rate of both modes i = 1,2 increases 
with increasing B held regardless of the sign of a;*, and thus the total production rate also 
monotonically increases. Furthermore, we again obtain the angle dependence symmetric 
with respect to 633 = n 12 in the production rate. 

Next we consider the case with |eE| 7 ^ 0 and \eB\ = 0. In this case, and Fjj ■ Fjj 


become [see Eq. (33) 


F-j = 2 [-{eQqjfE"^ - 2guieQq,£Ecos6sE\ , 
IFij-Fjj = -4:[{gUif£B + gUieQq^BEcosOgE] . 


(52) 
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FIG. 5: Light (up) quark production rate in a S field as a function of which stands for 

9qb with a parallel configuration of £ and B. We take strengths of chromo-electromagnetic fields 
as gB = g£ = 1 GeV^. 

In this expression, we have used cos9be = cosOse- Again, the results are the same as 
those of the case where B is antiparallel to £. Figure 6 shows the quark production rate 
as a function of the angle 9se with several values of strength of the electric held. As 
the electric held increases, the production rate increases for whole angle regions. This 
can be understood in a similar way to the previous case as follows: At 9£e = 7r/2, the 
factor aijbij = \£ij ■ Bij\ = {gUiYSB is independent of the electric held. As for each 
factor, Uij = ^\JA{eQqYE^ + IQ^gojiYEB — 2 [eQqYE‘^ decreases with increasing electric 
held, while bjj = + lQ{guJi)‘^£B + 2{eQg)‘^E^ increases. These behaviors 

are opposite to those of the previous case with \eE\ = 0 and \eB\ ^ 0, and thus the 
production rate a.t 9 = 7r/2 monotonically increases. At 9£e = 0, tt, the quark production 


rate (43) can be rewritten as 

2,^Tfl /Iquark 


1 


Nc 






“ 1 , flgOOilB 

COth ' A, ' UTT 


—e 
n 




(53) 

1=1 n=l ^ 

where the strength of the combined electric held has been dehned as £ij = \gu}i£ + eQqE\ 
for 9£e = 0 and £ij = \gui£ — eQqE\ for 9£e = tt. In the case of \gu)i£\ > \eQqE\, when 
the chromo-electric held and the electric held are (anti)parallel to each other, 9£e = 0 
{9£e = tt), with cjj > 0 {oji < 0), the strength of the combined electric held £ij linearly 
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FIG. 6: Light (up) quark production rate in an E field as a function of which stands for 

9£e with a parallel configuration of £ and B. We take strengths of chromo-electromagnetic fields 
as gB = g£ = 1 GeV^. 


increases with increasing E field, and thus coth ( ^^riTr) monotonically increases. When 




^eE = 0 {9sE = vr) with cuj < 0 {ui > 0), the held strength Sij linearly decreases with 
increasing E held, and coth slightly decreases and approaches unity. Then, after 

summing over all the modes z, the production rate (53) at Ose = 0 {Ose = tt) monotonically 
increases with increasing E held. On the other hand, in the case of \guJiS\ < \eQqE\, the 
production rate of both modes i = 1,2 increases with increasing E held regardless of the sign 
of uji, and thus the total production rate also monotonically increases. From these results, 
we expect that strong EM helds created in the early stage of relativistic heavy-ion collisions 
would largely ahect quark productions from a glasma (chromo-EM helds) depending on the 
held conhgurations, and would thus possibly inhuence the formation of QGP. 


B. Weiss potential with electromagnetic fields 

In this subsection, we will investigate the ehects of EM helds on the conhnement- 
deconhnement phase transition by using the ehective potential of the Polyakov loop in 
the presence of EM helds. 

Prior to going into the details, let us briehy explain the ehective potential without external 
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fields being imposed. The one-loop calcnlation at finite temperatnre in SU(2) gange theory 
and in the massless fermion limit yields the effective potential for the temporal component 


of the gange field {C = |4l) as 


yWetop] ^ l/We,ss|q ^ I/S'IC] , 


^Weiss 


Weiss r 


where the YM and qnark parts are given, respectively, by 

vS‘“|C] = + 

CtfPl = -Y^r‘ + h=r‘‘c^(2-c=). 


(54) 

( 55 ) 

(56) 


This resnlt is called the Weiss potential. In Fig. 7, we show the Weiss potential 
and its breakdown. We see that in the YM part, the minima appear at C = 0 and C = 1, 
reflecting the center symmetry C* —)■ C* -|- 1 in SU(2). Thus, selecting one of the two minima 
spontaneously breaks the center symmetry. Since the system should be in the deconhned 
phase in the high-temperature region where a perturbative approach becomes valid, this 
result seems to be natural. The quark part of the effective potential explicitly breaks the 
center symmetry, and C = 0 and C = 1 are no longer degenerated. In the presence of the 
quark part, C = 0 is favored, which corresponds to the deconhned phase. We are now going 
to investigate how this picture is modihed by the presence of external EM helds. 


Now we come back to our most general results (26), (27), and (40). Taking the vanishing 


limit of the chromo-EM helds, T, S —)■ 0, but keeping the Polyakov loop and EM helds 
nonzero in the results , we obtain the ehective potential 


Ves[A4:,E,B] — 


J dx^ 
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Y 1 /-OO J °° 

E / ^{4-2}2E^' 
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EE 

i=l f=l■ 
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—e *™'‘'J'^(a/s)(b/s)cot(a/s)coth(b/s) 


'0 


X 


2^(-l)Vs^ 


AhfisW I 

t ’2 '4'^ > cos ———n 


n=l 


T 


(57) 


where 0 / and bf are just given by the EM helds as 
1 


0/ = 


Ff + {F,-F,r + Ff, bf 


b,= - 


Fj + ( F, . F,Y - F] , 


(58) 
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FIG. 7: Weiss potential as a function of C. Constant terms which are independent of C are 
subtracted. 

with Fj = 2{eQqjY{B‘^ — E"^) and Fj ■ Ff = —A^eQq^YE ■ B. The factor f)/(s) is given by 

^2 _j_ g2 

‘^/(^) = S .U 2 Q/cot(a/g) + ^2 , J b/coth(b/g), (59) 

where Cj = {uaFj^){u 0 Fj^) = {eQq^YE'^ with = (1,0, 0,0). Here we have subtracted 
divergences appearing in the zero-temperature contribution, which are independent of A 4 ,. 


1 . ITeiss potential in magnetic fields 


Consider a pure magnetic held case, E ^ 0, B 0. Then, the effective potential reads. 


VeslAi, B] — 
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4 T 2 s cos 


h=l 

POO 


n=l 


gvhAi 


T 


-n 


j poo 1 2 

/ -^^~''"^^\e\Qq^\Bs)coi{e\Qq^\Bs) 

i=l /=! "'0 
00 

2E(-1)’ 


x; 


cos 


n=l 


gujiA4 
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-n 


(60) 


We rewrite the proper time integrals in two steps. Recall that the integral should be dehned 
with an inhnitesimally small number 6 which makes the contour slightly inclined to avoid 
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the poles along the real axis (in the second term). Then we can easily change the contour 
from [0, cxd] along the real axis to [— zcxd,0] along the imaginary axis (the Wick rotation), 
since there is no pole along the imaginary axis. Finally, by renaming the variable s as —ia, 
we obtain the following representation with integrals dehned by real functions®: 

Wff[A,5] = I 3t^cos( 


h=l 
Nc Nf 


(T"’ 


n=l 
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da _m 2 




'^‘^f'"ie\Qq^\B)coth{e\Qq^\Ba) 


i=i f=i 
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”e 4T2a COS 


n=l 


guiAi 


-n 


(61) 


For simplicity, we shall restrict ourselves to W = 2, which provides us with all the essential 
features of the perturbative effective potential in the presence of EM helds. In this case, the 
eigenvalues cuj and Vh are simply given by Ui = ±1/2 and Vh = 0, ±1. The effective potential 
reads. 


V,s[C,B\ = -^ir'=T*+?7rW(l-C)'^ 


(62) 


Nf 
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da _m 2 


27^2 i —/ j ijz 

f=l 


(e\Qqj,\B)coth.{e\Qq^\Ba) ^(-l)”e it^cos (Cvrn) 


n=l 


The hrst line does not depend on the magnetic held and corresponds to the YM part Vym- 


This is nothing but the Weiss potential (55) 05]. The second line corresponds to the quark 
part Yquark, and the integral and summation over n can be easily performed numerically. 
From now on, we further restrict ourselves to the one havor / = 1 with the electric charge 
Qqj = 1 for simplicity. Now, analytic expressions are available in two limiting cases: One is 
the i? —?• 0 and niq —?■ 0 limit, where the quark part of the effective potential is reduced to 


that of the Weiss potential (55): 


+ tx=r‘‘C"(2 - C=) = qsj [C] , (63) 

The other is the strong magnetic held limit: eB S> where the quark part can be written 


as 


EquarkiO, S] — —2 


[eB) 


TT^ 


vr 


— 
12 


(Ott) 


( 64 ) 


The second line of Eq. (611 coincides with Eq. (B.6) in the appendix of Ref. [T^ . 
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FIG. 8: Quark part of the effective potential as a function of C for several values of magnetic fields. 
X and y are given as x = m^/T^ and y = eB/T"^, respectively. 


Figure 8 shows the magnetic held dependence of the quark part of the effective potential 


which is given by the second line of Eq. (62). Here, we show only one havor contribution 


with X = = 0.5. An important observation is that as the magnetic held increases, the 

explicit breaking of the center symmetry is enhanced, and C = 0 (deconhned phase) becomes 
more stable. This is qualitatively consistent with the analytic representation at strong 


magnetic helds [see Eq. (64)] in that the potential value at G = 0 becomes more negative and 
the rising behavior becomes steeper with increasing magnetic held. The enhancement of the 
center symmetry-breaking ehects due to increasing magnetic held indicates that the quark 
loop interacting with magnetic helds can be one of the important sources for reducing the 
(pseudo)critical temperature Tc of conhnement-deconhnement phase transition, as observed 
in recent lattice QCD simulations [19]. In the last part of this subsection, we will see within 
a phenomenological model that this is indeed the case. 
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2. Weiss potential in electric fields 


In the case of a pure electric field, -B —)■ 0 and E 0, the situation is a bit subtle. The 
effective potential of the quark part can be written as 

Nf 

1 ^ ds ■ 2 

Vqni,rk[AA, E] = ^ (c||Es) coth (c||Es) 


/=1 
oo 


n=l 




(65) 


Note that we cannot reach this result from Eq. (60) by replacing B with iE, un¬ 
like the zero-temperature contribution. This is due to the form of the factor f)/(s) = 
(e|Qg^|E's)coth(e|Qg^|Es) in the exponential. Because of this factor, the full calculation 
(even numerical evaluation) is rather difficult. Furthermore, since there are singularities 
(poles) on the imaginary axis, we cannot perform the Wick rotation of the proper time s, 
unlike the Weiss potential in magnetic helds. To avoid these difficulties, we expand the 
effective potential with respect to the electric held. Using a:cothx ~ 1 -f x^/3 ■ • •, we get 
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( 66 ) 


At this stage, we can perform the Wick rotation for the proper time s. Then, the effective 
potential reads 

1 .00 
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The systematic expansion with respect to the E held is possible, and the integral and sum 
can be performed numerically at each order. 

In Fig. 9 we show the electric held dependence of the quark part of the ehective potential. 
From this hgure, we see that the electric held decreases the explicit breaking of the center 











FIG. 9: Quark part of the effective potential as a function of C for several values of electric fields. 
X and y are given as x = m^/T^ and y = eE/T"^, respectively. 

symmetry. This is completely opposite to the B dependence of the effective potential. Thus, 
we expect that Tc increases with increasing E field and approaches the Tc of the pure YM 
theory. 


3. Phenomenological analysis on Tc{B) 

We have seen that imposing magnetic fields enhances the explicit breaking of the center 
symmetry. What we have evaluated is a perturbative contribution (in the sense that we 
assume that the coupling is small enough), and thus we discussed how the Weiss potential 
(that is also evaluated in a perturbative framework) is modified in the presence of the EM 
fields. Within this perturbative calculation, we are not able to approach the region where 
phase transition will take place. Indeed, even if the quark part of the effective potential 
depends on the magnetic fields V)juark[C', 5], the total effective potential Ves\C, B] = Vym[C] + 
Vquark[C', .B] selects the center broken state (7 = 0, and thus confinement-deconhnement 
phase transition never occurs within this perturbative framework. However, recall that the 
magnetic field can affect the effective potential of the Polyakov loop only through the quark 
loop at leading order. Therefore, we expect that even the perturbative evaluation of the 
quark part V)juark[(7, 5] can make sense if combined with some nonperturbative effective 
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potential [C] for study of the effects of magnetic fields on the phase transition. Here 

we discuss whether this is indeed the case. 

Let us introduce a simple model of a gluonic potential reproducing conhnement- 
deconfinement phase transition, 

U[C] = -^a(T)$2 + 6 (T) In [l-6<h2 + 8<h3-3$^] ( 68 ) 

with 


a(T) = o„ + o,(ro/r) + a2{TolTf, b(T) = h(T„/T)' 


(69) 


Now, we consider the Nc = 3 case. Here the parameters are oq = 3.51, ai = —2.47, 02 = 
15.2, 63 = —1.75, and Tq = 270 MeV, which are hxed to reproduce the quenched lattice QCD 
results [55]. Instead of VVm, we employ this phenomenological potential ( 68 ) and combine 


it with V)juark[C',-B]. In this way, we can study how the temperature dependence of the 
Polyakov loop changes with magnetic helds. Notice that the quark part of the perturbative 
effective potential V)juark[C', B] with iVc = 3 is the same as that of the one with Nc = 2, since 
the quark with 0)3 = 0 does not contribute to the potential. Therefore, we can use the same 


potential evaluated in the second line of Eq. (62). The result is shown in Fig. 10. In this 


analysis, we have used Ui = ±1/2, 0 and a constituent quark mass rUg = 350 MeV. Thanks 
to the explicit center symmetry breaking, the Polyakov loop increases with increasing B 
field, in particular below the phase transition temperature, which eventually brings about 
decreasing pseudocritical temperature Tc{B) < Tc{B = 0). This result is very encouraging, 
but obviously we need to couple quark dynamics to the gluon dynamics to understand the 
effects of magnetic helds on the actual phase transition. 

Very recently, the inverse magnetic catalysis of the chiral sector, namely the decrease of 
the critical temperature of the chiral phase transition, has been reproduced from functional 
approaches including the Dyson-Schwinger equations and the functional renormalization 
group [241 [25] . Once the inverse magnetic catalysis of the chiral sector occurs, dynami¬ 
cal quark masses decrease with increasing magnetic held around T^. Then, the quark loop 
contribution is enhanced, and thus the ehect of the explicit center symmetry breaking be¬ 
comes larger. Therefore, the inverse magnetic catalysis of chiral sector would support the 
decreasing of the Tc of conhnement-deconhnement phase transition through the quark loop. 
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FIG. 10: Temperature dependence of the Polyakov loop for different values of magnetic field. 

IV. SUMMARY AND CONCLUSION 

In the present paper, we analytically derived the Euler-Heisenberg action for QCD+QED 
in the presence of the Polyakov loop, called the Euler-Heisenberg-Weiss action, by using the 
Schwinger proper time method. The effective action contains EM helds and chromo-EM 
helds as well as the Polyakov loop in a nonlinear form and reproduces the known one-loop 
effective actions for QED, QCD, QCD-I-QED, and also the Weiss potential for the Polyakov 
loop in appropriate limits. 

As an application of our effective action, we investigated quark pair productions under 
strong EM helds and chromo-EM helds. Using the ehective action of the quark part at zero 
temperature, we derived the formula describing the quark pair production rate in arbitrary 
conhgurations of the QCD helds and QED helds. In particular conhgurations, EM helds 
enhance the quark pair productions induced by chromo-EM helds. This indicates that strong 
EM helds created in relativistic heavy-ion collisions would largely ahect quark-antiquark 
pair productions from a glasma and thus could give sizable contributions to the formation 
of QGP. 

We also studied the perturbative ehective action of the Polyakov loop in the presence 
of strong EM helds. We found that the magnetic (electric) held enhances (reduces) the 
explicit center symmetry breaking through the quark loop. This indicates that the Polyakov 
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loop increases as the magnetic field increases, and thus the (pseudo)critical temperature 
of confinement-deconfinement phase transition decreases. In contrast, the electric field 
would raise the critical temperature. In order to demonstrate this, we combined the quark 
part of our perturbative effective potential with a simple model which can reproduce the 
confinement-deconfinement phase transition. The resultant Polyakov loop indeed increases 
with increasing B field, and then (pseudo)critical temperature decreases. This result is 
consistent with recent lattice data. Very recently, G. Endrodi investigated QCD phase 
transitions in unprecedentedly strong magnetic fields from lattice simulations of 1 -|- 1 -|- 
1-fiavor QCD He found strong evidence for a first-order confinement-deconfinement 

phase transition in the asymptotically strong magnetic field regions. In order to understand 
these lattice data, further nonperturbative analyses will be necessary. As a future work, we 
will extend the present work to nonperturbative analyses in terms of functional approaches. 
The inclusion of the chiral sector (quark-quark interaction mediated by gluons) will also be 
an important ingredient in the future work. 
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